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ABSTRACT: 

Neuroblastoma, the most common extracranial solid tumor of childhood, arises 
from neural crest cell precursors that fail to differentiate. Inducing cell differentiation 
is an important therapeutic strategy for neuroblastoma. We developed a direct 
functional high-content screen to identify differentiation-inducing microRNAs, in 
order to develop microRNA-based differentiation therapy for neuroblastoma. We 
discovered novel microRNAs, and more strikingly, three microRNA seed families that 
induce neuroblastoma cell differentiation. In addition, we showed that microRNA seed 
families were overrepresented in the identified group of fourteen differentiation- 
inducing microRNAs, suggesting that microRNA seed families are functionally more 
important in neuroblastoma differentiation than microRNAs with unique sequences. 
We further investigated the differentiation-inducing function of the microRNA-506-3p/ 
microRNA-124-3p seed family, which was the most potent inducer of differentiation. 
We showed that the differentiation-inducing function of microRNA-506-3p/microRNA- 
124-3p is mediated, at least partially, by down-regulating expression of their targets 
CDK4 and STAT3. We further showed that expression of miR-506-3p, but not miR- 
124-3p, is dramatically upregulated in differentiated neuroblastoma cells, suggesting 
the important role of endogenous miR-506-3p in differentiation and tumorigenesis. 
Overall, our functional screen on microRNAs provided the first comprehensive analysis 
on the involvements of microRNA species in neuroblastoma cell differentiation and 
identified novel differentiation-inducing microRNAs. Further investigations are 
certainly warranted to fully characterize the function of the identified microRNAs in 
order to eventually benefit neuroblastoma therapy. 



INTRODUCTION 

Neuroblastoma is the most common solid tumor of 
infancy and the most common extracranial solid tumor 
of childhood, accounting for more than 7% of childhood 
cancers and 15% of cancer-related childhood deaths [1,2]. 
Neuroblastoma arises from the neural crest cell precursors 
of the sympathetic nervous system that fail to differentiate 



[2, 3] - this provides the basis for differentiation therapy, 
an approach to induce malignant cells to differentiate 
into mature cells, thereby leading to cell growth arrest 
and apoptosis [2, 4-6]. However, only a limited number 
of differentiation agents have been successfully used to 
treat neuroblastoma. The differentiation agent 13-cis- 
retinoic acid (RA) is currently the standard of care 
for post-remission maintenance therapy in high-risk 
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neuroblastoma [2]. Although such treatment has resulted 
in a significant increase in patient survival, more than 
50% of the treated patients still develop recurrence [7, 
8]. Such poor outcomes demand the development of new 
differentiation agents. Unfortunately, the mechanisms that 
result in the loss of differentiation ability of neuroblastoma 
cells are poorly understood, which poses an obstacle to 
such development. Therefore, identifying additional 
differentiation agents largely relies on the discovery of 
new targetable biological molecules that play critical roles 
in neuroblastoma differentiation. 

High-throughput screening approaches significantly 
facilitate the discovery of novel anti-cancer drugs 
and drug targets. More recently, high-content screens 
(HCSs) based on automated cell imaging have been 
developed. However, current HCSs generally either 
use genetic engineered cell lines expressing fluorescent 
signals or involve staining of fixed cells [9, 10], which 
are generally time-consuming and consequently limit 
their broad applications to drug discoveries. We aim to 
develop a non-fluorescence, live-cell based HCS approach 
for identifying neuroblastoma differentiation- inducing 
agents. Neurite outgrowth, which is easily detectable 
under the microscope, is a well-recognized morphological 
differentiation marker of neuroblastoma cells in vitro [11- 
14]. While undifferentiated cells usually show no visible 
neurites, fully differentiated neuroblastoma cells form 
neurites that are four to five times the length of the cell 
body. This differentiation trait facilitates the design of a 
functional HCS assay to directly identify substances that 
induce neuroblastoma cell differentiation. In this study, we 
developed such an approach and applied it in a screen for 
microRNAs (miRNAs) that induce differentiation. 

miRNAs are endogenously expressed small RNAs 
that play a critical role in tumorigenesis [15-19]. The 
therapeutic potential of either exogenously increasing 
cellular miRNAs levels with synthetic miRNA mimics, or 
inactivating endogenous miRNAs with synthetic miRNA 
inhibitors has been demonstrated in previous studies [20- 
22]. The role of miRNAs in neuroblastoma differentiation 
and tumorigenesis has been implicated[23-31], which 
suggests the potential of developing novel miRNA- 
targeting approaches to neuroblastoma differentiation 
therapy [32], and warrants a comprehensive understanding 
of the involvement of miRNAs in neuroblastoma cell 
differentiation. However, there has been no concerted 
effort to comprehensively investigate the functions of 
the miRNA species in neuroblastoma differentiation. By 
applying the HCS that we developed, we investigated 
the recently identified human miRNAs and identified 
differentiation-inducing miRNAs that have not been 
discovered previously. 



RESULTS 

A HCS approach for measuring neuroblastoma 
cell differentiation is developed based on neurite 
quantification. 

Neurite outgrowth is well recognized as a 
morphological hallmark of neuroblastoma cell 
differentiation in vitro [11-14]. This facilitates 
the development of a HCS approach to identify 
differentiation-inducing agents based on quantification 
of neurite outgrowth. Neuroblastoma cell line BE(2)-C 
shows easily detectable neurite outgrowth upon induced 
differentiation by all-trans retinoic acid (ATRA). As 
shown in Figure 1 A, ATRA (b) induces dramatic neurite 
outgrowth in BE(2)-C compared to control (a), and the 
neurites and cell body area can be clearly defined (c, d). 
Quantification (Figure IB) shows that ATRA significantly 
increases the relative neurite length compared to control. 
In addition, ATRA induces neurite elongation in both time- 
and dose-dependent manners (Figure IC-D and Suppl. 
Figure 1). Correspondingly, ATRA decreases cell viability 
(Figure IE), stimulates expression of neuroblastoma 
differentiation markers (i.e., growth associated protein 43 
(GAP43), neuron specific enolase (NSE) and P-TUBULIN 
III) [33-35], inhibits expression of cell proliferation 
markers (i.e., PCNA and Ki67), and increases expression 
of apoptosis markers (i.e., cleaved CASPASE 3 and 
PARP) (Figure IF) in dose-dependent manners. These 
results indicate that neurite length is a reliable quantitative 
marker of BE(2)-C cell differentiation, and therefore can 
be used to compare the efficacy of differentiation-inducing 
agents. This was the basis of our HCS protocol (Suppl. 
Figure 2) for identifying novel differentiation-inducing 
miRNAs. 

HCS identifies novel miRNAs that induce 
neuroblastoma cell differentiation. 

Using our HCS protocol, we screened a library of 
miRNA mimics (Dharmacon) in BE(2)-C cells. Figure 
2A shows the neurite length distribution (grey histogram) 
associated with individual miRNA mimics. Replicate 
screens for one library plate from two independent 
transfections show that the results are highly reproducible 
(Figure 2B) (R=0.95,;?<0.0001), supporting the reliability 
of the screen. As shown in Figure 2A, the neurite length 
distribution following treatment with miRNA mimics is 
asymmetric; a small number of miRNAs are identified as 
dramatically increasing neurite lengths on the far right side 
of the distribution. In order to examine the neurite length 
distribution of unaffected cells, 13 plates of untreated 
BE(2)-C cells were measured (Figure 2C). Kolmogorov- 
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Figure 1: Neurite length is a quantifiable differentiation marker of BE(2)-C cells. 2,500 cells were plated in 96-well plates 
and cultured overnight. Cells were then treated with ATRA or carrier (DMSO, control) and placed into the IncuCyte for detecting neurite 
outgrowth. 9 images were taken from each well to allow for statistical analysis. Relative neurite length is defined as neurite length per 
cell body area. A, ATRA induces neurite outgrowth. Shown are representative phase-contrast images for cells treated with (a) carrier or 
(b) ATRA for 5 days, and (c, d) the same images analyzed to define neurites (pink) and cell body areas (yellow). B, Quantification shows 
that ATRA significantly increases the relative neurite length compared to control. C, Relative neurite lengths increase in a time-dependent 
manner during ATRA-induced cell differentiation. Neurite lengths were normalized to the starting time point (0 h). D, Dose-dependent 
effect of ATRA on neurite outgrowth. Shown are the results after treating with ATRA for 5 days. E, Dose-dependent effect of ATRA on 
cell viability. Cells were treated with different concentrations of ATRA, and cell viability was determined after 5 days. F, Dose-dependent 
effect of ATRA on the protein expression levels of cell differentiation markers GAP43, NSE and P-TUBULIN III, cell proliferation markers 
Ki67 and PCNA, and apoptotic markers cleaved CASPASE 3 and PARP, with GAPDH protein levels used as a loading control. Cells were 
treated with ATRA as above, and protein levels were determined by Western blots after 5 days. **,/?<0.01; ***,/?<0.001; ****,/?<0.0001. 
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Figure 2: HCS screening of miRNA mimics identifies 
miRNAs that induce neurite outgrowth in BE(2)-C cells. 

A-B, 2,500 cells were reverse-transfected with 25 nM miRNA 
mimics in 96-well plates. After 4 days of transfection, relative 
neurite lengths were quantified as above. A, The distribution of 
post-normalized neurite length measurements (grey histogram) 
for individual miRNA mimics from the screen and the fitted 
Gamma distribution curve (red line, a=3.27, P=0.444) generated 
by the method described in Suppl. Figure 3. B, Correlation 
between two independent neurite length measurements for one 
library plate. Cells were transfected with 25 nM mimics in two 
independent 96-well plates. Normalized neurite lengths were 
measured as above. Correlation of neurite length between the 
two plates was analyzed using two-tailed Pearson Correlation 
with p<0.05 considered significant. C, The distribution of the 
neurite lengths for untreated BE(2)-C cells. 2,500 Cells were 
plated into each well in 96 well plates. After 4 days, neurite 
lengths were analyzed as above. Shown are the neurite length 
distribution histogram (grey), the empirical density curve (black 
line), and the fitted Gamma distribution curve (red line, a=30.23, 
P=0.033). Using the P<0.05 threshold, Kolmogorov-Smimov 
goodness-of-fit test for Gamma model validity did not reject the 
null hypothesis (p=0.16), which support that the neurite length 
distribution fits Gamma model. 
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Figure 3: Characterization of the effect of the top 5 neurite-inducing miRNA mimics on cell differentiation and growth 
in multiple neuroblastoma cell lines. A-B, Effects of the identified top 5 miRNA mimics on neurite outgrowth. BE(2)-C cells were 
transfected with 25 nM miRNA mimics or control for 4 days, and relative neurite lengths were quantified as above. Shown are representative 
cell images analyzed to define neurite and cell body areas (A) and neurite length quantifications (B). C, Effects of the 5 miRNA mimics 
on expression of differentiation markers, with CALNEXIN protein levels used as a loading control. BE(2)-C cells were transfected with 
25 nM of the indicated miRNA mimics or control, and proteins levels were examined after 4 days. D, Effects of the 5 miRNA mimics on 
cell proliferation rate. E, Effect of the 5 miRNAs on neurite outgrowth in multiple cell lines. Cells were transfected with 25 nM miRNA 
mimics or control, and neurite lengths were quantified as above after 5 days. F-G, Effects of miR-506-3p (F) and miR-124-3p (G) mimics 
on expression of differentiation markers in multiple cell lines. *,/?<0.05; **,/?<0.01; ***,/?<0.001; ****,/7<0.0001. 
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Table 1: Fourteen miRNA mimics identified from HCS as inducing neurite outgrowth using a FDR threshold <0.01. Shown are the 
(a) miRNA name, (b) fold change of neurite length relative to unaffected cells, (c) p value, (d) FDR, (e) mature sequences of the miRNAs 
with seed sequences underlined and seed family grouping indicated in parentheses, and (f) references that have reported the involvements 
of the corresponding miRNAs in differentiation. 



miRNA 


(b) Neurite length 
(MeaniSD) 




\ FDR 






hsa-miR-124-3p 


14.55±4.75 


l.llE-11 


1.37E-08 


UAAGGCACGCGGUGAAUGCC (1) 


27,37 


lloCl llllXV Iw/w^U 


1 4 02±4 29 


3 26E-11 


2 02F-0R 


T J AT IGGCT JI JI JI ir AT JI JCCT J AT JGI JGA 




nsa-miK-jUo-Jp 


1 A A/|_|_0 A 1 


l.Ujn-U / 


/I l/IU 


T T A A r^r^r^ a /^/^/^t tt t/^t t/^ a r^T t a a ^^\ 




hsa-miR-34a-5p 




o 1 1 T7 f\n 


Q Q 1 AC 


ULrLrCALrULrUCU UALrCULrLrU ULrU (z) 




nbd-iniiv- 1 ujd- jp 




J. JOJJ/-U / 


O.O 1 JJ/-UJ 


A c^c A nr^ A T TT THT T A A nnnr^T T A T in A 


zo 


Jisa-miK-4jUD-Jp 


y.zo±z.o4 


/I q/ith c\n 
4.y4r/-U / 


1 noT^ n/i 
1 .Uzr/-U4 


T TT innn a t t/^ a t tt tt tt inr^ a t ir^r^ a t t a 






Q 1 ^+9 ZL7 


^ O^Th 07 
O.UJJ_/-U / 


1 n7F OZL 






hsa-miR-2110 


9.00±2.40 


8.21E-07 


1.27E-04 


UUGGGGAAACGGCCGCUGAGUG 




hsa-miR-34b-5p 


8.4U2.16 


2.66E-06 


3.66E-04 


UAGGCAGUGUCAUUAGCUGAUUG 




hsa-miR-107 


7.60±1.98 


1.27E-05 


1.55 E-03 


AGCAGCAUUGUACAGGGCUAUCA (3) 


36 


hsa-miR-3714 


7.56±1.91 


1.37E-05 


1.55E-03 


GAAGGCAGCAGUGCUCCCCUGU 




hsa-miR-449b-5p 


7.40±1.85 


1.89E-05 


1.95E-03 


AGGCAGUGUAUUGUUAGCUGGC (2) 




hsa-miR-137 


7.03±1.36 


3.83E-05 


3.65E-03 


UUAUUGCUUAAGAAUACGCGUAG 


19, 34 


hsa-miR-3937 


6.93±1.05 


4.61E-05 


4.08E-03 


ACAGGCGGCUGUAGCAAUGGGGG 





Smimov goodness-of-fit test for Gamma model validity 
indicates that the neurite length distribution fit Gamma 
model (p=0.16), which informs us to use this model 
(red line in Figure 2A) to assess the effect of individual 
miRNAs, as illustrated in Suppl. Figure 3. Fourteen 
miRNA mimics were identified as significantly increasing 
neurite length using False Discovery Rate (FDR) threshold 
<0.01 (Table 1). Using the same threshold, 0 hits were 
generated from untreated cells, indicating the specificity 
of the analysis approach. Among the 14 miRNAs, several 
were related to neuroblastoma cell differentiation in 
previous studies (Table 1) [27, 28, 36, 37], demonstrating 
the sensitivity of our HCS approach. 

We next focused on characterizing the 
differentiation-inducing function of the top 5 microRNAs 
that are most potent in inducing neurite outgrowth (Figure 
3A-B). We showed that, comparing to control, the 5 
miRNAs induce expression of differentiation markers and 
dramatically decrease BE(2)-C cell growth rate (Figure 
3C-D), demonstrating that true cell differentiation and 
growth arrest are induced. We further investigated the 
5 miRNAs in additional neuroblastoma cell lines with 
different genetic backgrounds (Suppl. Table 1). Among the 
5 miRNAs, miR-506-3p and miR-124-3p, which belong 
to the same seed family (defined as a group of miRNAs 
that share common seed sequence. Table 1), have the most 
potent effect on neurite growth (Figure 3E and Suppl. 
Figure 4). Correspondingly, the two miRNAs dramatically 
induce expression of differentiation markers in all the 
tested cell lines (Figure 3F-G). We further show that 
miR-506-3p and miR-124-3p mimics significantly reduce 



the ability of BE(2)-C to form colonies, indicating their 
long-term capacity to inhibit cell proliferation (Figure 
4A-C). To exclude the possibility that the induced cell 
differentiation and growth arrest are caused by off-target 
effects of the specific chemical designs of miRNA mimics, 
we examined the effect of miR-506-3p and miR-124-3p 
precursors (Ambion) on differentiation. The miRNA 
precursors are partially double-stranded RNAs designed 
to mimic the functions of the endogenous miRNAs, in 
contrast to the fully complementary double- stranded 
design of miRNA mimics. Figure 4D-E shows that the 
precursors significantly induced neurite outgrowth, 
recapitulating the results with miRNA mimics. This 
indicates that the differentiation-inducing function of 
miR-506-3p and miR-124-3p mimics is unlikely caused 
by off-target effect. Overall, the above results demonstrate 
the general and potent effect of miR-506-3p/miR-124- 
3p on inducing differentiation, suggesting the potential 
of restoring miR-506-3p and miR-124-3p expressions 
as a novel differentiation therapeutic strategy to treat 
neuroblastoma. 

To further examine the potential pathophysiological 
relevance of endogenous miR-506-3p and miR-124-3p in 
regulating neuroblastoma differentiation, we measured 
their expression levels in differentiated BE(2)-C cells. 
Figure 4F shows that expressions of both miRNAs are 
significantly increased in differentiated cells. However, 
the overexpression of miR-506-3p (3 80.1 ±48. 4 fold) is 
much more dramatic than that of miR-124-3p (8.1±0.4 
fold). These results indicate that endogenous miR-506- 
3p expression in neuroblastoma cells is highly regulated 
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Figure 4: Characterization of the tumor suppressive function of miR-506-3p/miR-124-3p family. A-C, Colony formation 
assay as a function of miR-124-3p and miR-506-3p mimics. BE(2)-C cells were transfected with the 25 nM of the indicated oligos and 
colony formation was examined as above. Shown are (A) plate images of colony formation and quantified colony numbers (B) and sizes 
(C). D-E, Effect of miR-124-3p and miR-506-3p precursors on neurite outgrowth. BE(2)-C cells were transfected with or without 10 nM 
miRNA precursors. Shown are (D) Representative images analyzed to define neurites and cell body areas after 4 days transfections and 
(E) quantification of relative neurite lengths. F, Effect of cell differentiation on endogenous expression of miR-506-3p and miR-124-3p in 
BE(2)-C cells. G, Dose-dependent effect of miR-506-3p mimic on neurite outgrowth in BE(2)-C cells. Cells were transfected with different 
concentrations of miR-506-3p mimic, and relative neurite lengths were quantified as above after 4 days. Neurite lengths were normalized to 
control (0 nM). H, Dose-dependent effect of miR-506-3p mimic and cis-RA on cell viability in BE(2)-C cells. Cells were transfected with 
different concentrations of miR-506-3p mimic for 4 days or treated with different concentrations of cz^'-RA for 5 days, and cell viability was 
determined as above. *,;?<0.05; **,/><0.01; ***,;?<0.001; ****,/><0.0001. 
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and suggest that, relative to miR-124-3p, endogenous 
miR-506-3p is likely to be a more dominant driving force 
in controlling cell differentiation. We further observe 
that miR-506-3p mimic induces neurite outgrowth and 
decreases cell viability in a dose-dependent manner, and 
that the induction of neurite outgrowth is significant at 
a concentration as low as 0.1 nM (Figure 4G-H, Suppl. 
Figure 5). Figure 4H also shows that the cytotoxic effect 
of miR-506-3p mimic is much more potent (reducing cell 
viability to 6.3±0.3% when reaching plateau) than 13-cis 
retinoic acid (cis-RA) (42.6±3.7%), a differentiation- 
agent currently used to treat neuroblastoma patients [7]. 
Altogether, these results demonstrate the potent function 
of miR-506-3p in inducing differentiation and in reducing 
cell survival and growth. The role of miR-506-3p in 
neuroblastoma tumorigenesis has not been investigated 
previously; we are the first to discover its differentiation 
inducing function in neuroblastoma. 

From the screen, we also identified miRNAs 
that significantly reduce cell survival but do not 
induce differentiation (Suppl. Figure 6). This indicates 
that miRNAs modulate neuroblastoma cell survival 
and growth through distinct mechanisms, and cell 
differentiation is not a prerequisite for cell death or growth 
arrest, which supports the functional specificity of the 
identified differentiation-inducing miRNAs in regulating 
differentiation. 

miRNA seed families that are potent inducers of 
cell differentiation are discovered from HCS. 



miRNA seed sequences are central in determining 
their target genes [38-41]. We therefore first analyzed the 
seed sequences of the identified differentiation-inducing 
miRNAs. Surprisingly, as shown in Table 1, three miRNA 
seed families, which accounts for 7 miRNAs, are identified 
within the 14 miRNAs. Enrichment analysis by random 
permutation shows that the probability that >7 non-unique 
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Figure 5: Enrichment analysis of seed families in a set 
of 14 miRNAs by random permutation. Grey bar shows 
the probability (p=2.2XlO'^) that at least 7 non-unique seed 
sequences appear in a randomly selected set of 14 miRNAs. 



seed sequences appear in a set of 14 randomly drawn 
miRNAs from the miRNA mimic library is p=2.2XlO'^ 
(Grey bar in Figure 5), indicating that miRNA seed- 
sequences families are significantly overrepresented in 
the identified 14 miRNAs. Further investigation of the 
remaining 7 miRNAs shows that the seed sequences of 
miR-135b-5p (shares seed sequence with miR-135a-5p), 
miR-34b-5p (shares with miR-449c-5p and miR-2682- 
5p) and miR-450b-3p (shares with miR-769-3p) are not 
unique. In addition, seed sequence family 2 includes 
another miRNA, miR-34c-5p, which is not identified in 
the top 14 candidates. Close examination of the screen 
results indicates that miR-34c-5p, miR-135a-5p and miR- 
449c-5p also increase the neurite lengths (Suppl. Figure 
7), ranking as 17*, 43* and 128* in the screen (miR-2682- 
5p was not in the library), although the effect of miR- 
449c-5p did not reach statistical significance using p< 0.05 
threshold. miR-769-3p does not induce neurite outgrowth, 
and, however, is the only exception among the identified 
differentiation-inducing seed families. 

The identified differentiation-inducing miRNAs 
are predicted to target distinct spectra of genes 
involved in neuroblastoma differentiation. 



The above observations lead us to identify potential 
miRNA targets based on seed sequence matches. We 
first identified 48 genes that have been previously 
demonstrated to regulate neuroblastoma differentiation 
(Suppl. Table 2.) using Ingenuity Pathway Analysis (IPA) 
(Ingenuity System). Not surprisingly, further IPA miRNA 
target analysis indicates that each of the miRNAs/seed 
families is predicted to target multiple genes involved in 
neuroblastoma differentiation (Figure 6), among which the 
miR-506-3p/miR-124-3p family is predicted to target 10 
of the 48 genes. The results also show that, although the 
predicted targets of these miRN As/seed-sequence families 
overlap, each miRNA/seed-sequence family has a unique 
targetome, suggesting that they are likely to induce cell 
differentiation through distinct, but overlapping pathways. 

CDK4 and STAT3 play a role in mediating the 
differentiation-inducing function of miR-506-3p/ 
miR-124-3p family. 

We further investigated the targets of miR-506-3p/ 
miR-124-3p family. Table 2 shows that the expression 
changes of the 10 predicted target genes of miR-506-3p/ 
miR-124-3p induced by their overexpressions are almost 
identical (Table 2), further demonstrating that the seed- 
sequence is dominant in determining the miRNA function. 
The two miRNAs dramatically down-regulate two of the 
ten predicted targets, CDK4 and STATS. Figure 7 A shows 
the interactions of the two miRNAs with the predicted 
target sites in the 3'UTR of CDK4 and STAT3. We 
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validated their target sites in the 3'UTRs of CDK4 and 
STATS using luciferase reporter assays (Figure 7B-C). 
We further show that overexpressions of the two miRNAs 
down-regulate endogenous CDK4 and STATS protein 
levels (Figure 7D). Figure 7E-F shows that individual 
repression of CDK4 and STATS expression induces neurite 
outgrowth, and that their combined repression has an 
enhanced effect on neurite outgrowth relative to individual 
repression. These results indicate that CDK4 and STATS 
mediate, at least partially, the differentiation-inducing 
function of miR-506-Sp/miR-124-Sp, and suggest that the 
effect of miR-506-Sp/miR-124-Sp on cell differentiation 
are most likely mediated by concordantly down-regulating 
multiple target genes. 

DISCUSSION 

In this study, we developed a HCS approach to 
facilitate the discovery of novel differentiation-inducing 
agents for neuroblastoma. Several HCS approaches 
based on quantification of neurite outgrowth have been 
described [SS, 42-44]. However, quantifications of neurite 
outgrowth in these approaches were either based on cell 
lines engineered to express fluorescent reporters or involve 



Table 2: Changes of expression for the 10 predicted 
targets of the miR-506-3p/miR-124-3p family induced 
by miR-506-3p and miR-124-3p overexpression. 

Shown are the gene name and the fold change of mRNA 
expression levels (miRNA mimic vs. control) induced after 24 
h transfection with (1) miR-506-3p mimic and (2) miR-124-3p 
mimic. *, Shown are average expressions detected by multiple 
probes in the microarray. 



Gene 


Gene expression* 


(1) miR-506-Sp 


(2) miR-124-Sp 


CDK4 


0.S9 


0.4S 


STATS 


0.75 


0.75 


CEBPA 


0.88 


0.86 


PML 


1.01 


1.00 


RARG 


1.02 


l.OS 


BDNF 


1.04 


1.06 


THRB 


1.04 


1.06 


RHOQ 


1.07 


LIS 


SHCS 


1.08 


1.07 


AHR 


1.19 


1.17 



IGF1R 



RARG 




ESR1 



5GC 

(miR-50f3p/mlR^^3p) 

/ V^THRB 

•4 



AGGCAGU 
(miR-34b-5p) 



CDC42 
UAUUGCl>\ ^ 
{miR-137-3p) <i 



/ 

V 

CDKN1A 



AUGGCUU 
(miR-135b-5p) 

AAGGCAG 

(miR-3714) GGCAGUG 

(miR-34a-5p/miR-449b-5p/miR-449a) 



NFKBIA 



^ Enzyme 
^ Growth factor 
^ Kinase 

Nuclear Receptor 
CO Transcription Regulator 
^ Transmembrane Receptor 
O Unknown 
'iEi' microRNA 

— mlRNA target Relationship 
HI microRNA 
^1 Shared targets 
IH Unique targets 



Figure 6: The predicted differentiation-inducing targetome network for the identified 14 differentiation-inducing 
miRNAs. The 14 miRNAs were grouped into 1 0 seed-sequence groups. Predicted targets previously reported as involved in neuroblastoma 
differentiation were used to create the targetome network. 
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staining of fixed cells [33, 42-44]. These approaches are 
generally time-consuming. In addition, whether the extent 
of neurite outgrowth in neuroblastoma cells is a reliable 
marker to evaluate the potency of differentiation- inducing 
agents was not clearly characterized. Here we showed 
that, in neuroblastoma cell line BE(2)-C, neurite length of 
differentiated cells is quantifiable based on phase-contrast 
images of live cells and is a reliable marker of the extent 
of cell differentiation. This supports us to develop a HCS 
approach to examine neuroblastoma cell differentiation by 
quantifying neurite outgrowth, and we for the first time 
show that the neurite length distribution of untreated cells 
fit the Gamma model - this informs us to use this model 
to perform statistical analysis on neurite-based HCS. Our 
screen of a library of miRNA mimics demonstrates that 
our screening and statistical analysis methods are specific 
and sensitive for identifying differentiation-inducing 
agents. We expect that this approach will facilitate future 
discovery of novel differentiation- inducing drugs and drug 



targets for treating neuroblastoma. 

Our study provides the first comprehensive 
and direct functional analysis of miRNA species in 
inducing neuroblastoma cell differentiation. Our screen 
not only recapitulated several previous findings, but 
also identified miRNAs that have not been known to 
regulate neuroblastoma differentiation. For example, 
we are the first to identify the function of miR-506-3p 
in promoting neuroblastoma cell differentiation. Our 
further investigations show that endogenous miR-506-3p 
expression was dramatically increased in differentiated 
neuroblastoma cells. In addition, published data 
have indicated the tissue-specific expression of miR- 
506 in adrenal gland, the primary tissue of origin for 
neuroblastoma [45]. This evidence altogether implicated 
the important role of miR-506-3p in neuroblastoma 
pathogenesis. Overall, our comprehensive analysis of 
miRNAs in neuroblastoma cell differentiation is the first 
step towards elucidating the entire picture of miRNA 
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involvement in differentiation. We expect that further 
investigation on the identified individual miRNAs will 
provide a better understanding of the mechanisms of 
neuroblastoma tumorigenesis. 

In our HCS, several miRNAs that were previously 
reported to induce neuroblastoma cell differentiation were 
not identified as potent inducers of cell differentiation [29, 
30]. One of the possible explanations is the cell context- 
specific variation of the function of these miRNAs, in 
that our screen was conducted in a different cell line from 
those used in previous studies. Further studies are certainly 
warranted to compare the differentiation-inducing function 
of these previously reported miRNAs with the miRNAs 
identified in our study in a larger panel of neuroblastoma 
cell lines and in in vivo studies, and to compare their 
clinical relevance in neuroblastoma tumorigenesis in 
clinical studies. 

Another intriguing finding in our study is the 
overrepresentation of a set of miRNA seed families in 
our identified differentiation-inducing miRNAs. This 
observation is beyond providing convincing evidence 
supporting the reliability of our HCS approach. It leads 
us to an interesting speculation: the conservation of 
seed sequence among different miRNAs located at 
different genomic regions may be an important evolution 
trait that is selected by nature to guard the normal cell 
differentiation process during development; if one miRNA 
fails its expression in cells, expression of another seed- 
family member will perform similar function to prevent 
differentiation error. If this is true, we expect that miRNA 
seed sequence families are more likely to be involved in 
cell differentiation than miRNAs with unique sequences. 
There are certainly other possible explanations. Whether 
miRNA seed families are also overrepresented in other 
differentiation/development processes and what are the 
pathophysiological significances of this overrepresentation 
are certainly interesting questions to pursue in the future. 

Previous studies have demonstrated the promise 
of synthetic miRNA mimics as therapeutic agents in 
cancer treatment [20, 22]. Our identification of novel 
differentiation- inducing miRNA mimics therefore directly 
provides a group of novel candidates potentially applicable 
to treat neuroblastoma. Our identification of several seed- 
sequence families that are potent differentiation inducers 
reinforces the notion that miRNA seed sequences play 
the key role in defining their biological function, which 
provides the rationale for developing novel synthetic seed 
sequence-based oligos as differentiation-inducing agents. 
By replacing the nucleic acids in the non-seed positions 
while keeping the seed sequence unchanged, various seed 
sequence-based synthetic oligos can be designed and 
tested in vitro and in vivo, in order to identify the optimal 
design that has the most potent effect on neuroblastoma 
differentiation and has minimal non-specific cytotoxicity 
on normal cells and tissues. 

Currently, differentiation agents are limited to be 



used for post-remission maintenance therapy in high-risk 
neuroblastoma [2]. One of the reasons for this limitation is 
that the current available differentiation agents are not as 
potent as other anti-cancer agents in ablating cancer cells. 
However, our study show that miR-506-3p mimic reduced 
neuroblastoma cell viability to a much greater extent than 
cis-RA, raising the possibility of applying differentiation- 
inducing miRNA mimics to remission induction therapy. 
Further investigations in in vivo models are certainly 
warranted to examine the therapeutic efficacy of miR-506- 
3p mimics. 

In this study, we preliminarily analyzed the predicted 
targetomes of the identified top 14 differentiation-inducing 
miRNAs using informatics tools, and took a small step 
toward validating the predicted targets of miR-506-3p/ 
miR-124-3p family. Our investigation suggests that the 
differentiation-inducing function of a miRNA is likely 
mediated by concordantly down-regulating multiple 
targets. Our current analysis is limited to genes that are 
known to regulate neuroblastoma differentiation. It is 
highly likely that there are undiscovered targets of these 
miRNAs that play important roles in mediating their 
differentiation-inducing functions. Indeed, our expression 
array analyses indicate that miR-506-3p/miR-124-3p 
overexpression down-regulates many more of their 
predicted targets than we investigated in this study (data 
not shown). Whether these targets play roles in mediating 
their differentiation-inducing function warrants further 
investigation. 

In summary, we have established a HCS platform 
to screen for novel differentiation-inducing substances in 
neuroblastoma cells. We have identified novel miRNAs 
and miRNA seed families that induce neuroblastoma 
cell differentiation. Our study not only provides the first 
comprehensive understanding of the role of miRNAs in 
neuroblastoma differentiation, but also provides novel 
leads for developing miRNA-based differentiation agents 
for neuroblastoma treatment. We expect that further 
investigation based on our findings will not only promote 
better understanding of the mechanisms of neuroblastoma 
tumorigenesis, but may also provide new therapeutic 
strategies for treating neuroblastoma. 

MATERIALS AND METHODS 



Materials 

ATRA and cis-RA were from Sigma (St Luis, MO, 
USA). Dharmacon miRNA mimic library and individual 
miRNA mimics were from Thermo Fisher Scientific 
(Rockford, IL, USA). miRNA precursors were purchased 
from Ambion (Foster City, CA, USA). Rabbit anti-GAP43, 
anti-NSE, and anti-p-TUBULIN III were from Abeam 
(Cambridge, MA, USA). Rabbit anti-CALNEXIN, 
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anti-GAPDH and goat anti-rabbit secondary antibody 
conjugated with horseradish peroxidase (HRP) were 
from Santa Cruz (Dallas, TX, USA). Rabbit anti-PARP 
(cleaved), anti-CASPASE-3, anti-STAT3, and anti-CDK4, 
were from Cell Signaling (Danvers, MA, USA). Rabbit 
anti-Ki67 was from Millipore (Billerica, MA, USA). 

Cell lines. 

BE(2)-C cells were purchased from the ATCC. 
Other cell lines were obtained from the cell line repository 
at the Greehey Children's Cancer Research Institute. Cells 
were grown in DMEM/F12 supplemented with 10% fetal 
bovine serum. 

Detection and quantification of neurite outgrowth. 

Cells were plated and treated in 96-well plates. For 
detecting neurite outgrowth, cells were placed into ZOOM 
IncuCyte Imaging System (Essen Bioscience) and cell 
images were taken under 20X microscopic magnification. 
For detecting neurite outgrowth in a time- dependent 
manner, cell images were taken every 6 h. The neurite 
lengths associated with each treatment were calculated 
using the neurite definition defined for each specific cell 
line using the Neuro Track system (Essen Bioscience). 

Analysis of HCS data. 

The relative neurite length associated with cells in 
each well on the screen plates was determined as above. In 
order to allow for plate-to-plate comparison, neurite length 
associated with each well in each plate was first internally 
normalized to the mean of the corresponding plate, and 
multiple screen plates were then aggregated together to 
generate the neurite length distribution. The data was 
then further analyzed to determine the distribution of 
the unaffected cells as described below, and to identify 
differentiation-inducing miRNA mimics as illustrated in 
Suppl. Figure 3. 

Kolmogorov-Smirnov goodness-of-fit test for 
Gamma model validity. 

In order to examine whether the neurite length 
distribution of the untreated cells fit Gamma distribution, 
Gaussian kernel smoothing was first performed to 
generate the empirical density curve based on the neurite 
length distribution histogram of the untreated cells, and 
Gamma model parameters (a, P) were then estimated by 
minimizing the area in between empirical and Gamma 
distribution curves. Statistical significance of the fitness 
to Gamma distribution was examined by Kolmogorov- 
Smirnov goodness-of-fit test, with p<0.05 considered 



as rejecting the null hypothesis that the neurite length 
distribution fits Gamma distribution. 

Western blots. 

Cell lysates were prepared using RIPA buffer. 
Protein concentration was determined using the 
Pierce BCA assay (Thermo Fisher Scientific). For 
electrophoresis, equal amounts of cell lysate were resolved 
by SDS-PAGE and transferred to Immun-Blot PVDF 
membranes (Bio-Rad Laboratories). Membranes were 
blocked and probed with primary antibodies to specific 
proteins. Bound antibodies were detected with secondary 
antibodies conjugated with horseradish peroxidase (HRP) 
visualized by enhanced chemiluminescent (ECL) substrate 
(Pierce). 

Cell growth rate assay. 

Cells were plated in 96-well plates and were 
treated with specified conditions. Cells were placed into 
the IncuCyte imaging system and cell confluence was 
monitored every 6 h for 120 h. Cell growth curves were 
derived by comparing the cell confluences at different time 
points. 

Cell viability assay. 

Cells were plated in 96-well format and treated 
as specified. After culturing for 120 h, cell viability was 
determined using the CellTiter-Glo Luminescent Cell 
Viability Assay (Promega). 

Colony formation assay. 

Cells were transfected as specified and cultured 
overnight. 500 cells were re-plated in 10 cm dishes. After 
14 days, colonies were visualized by staining with 1% 
crystal violet. Colony numbers and sizes were analyzed 
using Image J (NIH, Bethesda, MD). 

Enrichment analysis of miRNA seed families. 

We use random permutation to examine whether 
miRNA seed families are significantly enriched in the 
identified set of differentiation-inducing miRNAs. The 
miRNA mimic library includes mimics for 1231 human 
miRNAs listed in miRBase 16.0. Among these miRNAs, 
900 miRNAs have unique seed sequences; 331 miRNAs 
share seed sequences with at least one other miRNA, 
constituting 125 seed sequence families. To calculate 
the probability of miRNAs from the same seed families 
randomly appears in a set of 14 miRNAs, 14 miRNAs 
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were randomly drawn from the 1231 miRNAs, and the 
number of non-unique seed sequences within the 14 
miRNAs is counted (X). 10^ permutations were run to 
determine the probability that at least X non-unique 
seed sequences appear in a randomly selected set of 14 
miRNAs (p(X)). 

miRNA target prediction and pathway analysis. 

miRNA target sites in 3'UTRs were predicted 
based on seed sequence complementarity and were 
identified using the IPA Pathway analysis function, which 
identifies any 7-nucleotide region (3 '-5') in a given 
3'UTR completely complementary to the seed sequence 
of a miRNA (2"^^ - 8* or - 7* nucleotides, 5 '-3') as 
a potential target site of this miRNA. To identify targets 
that potentially mediate the differentiation-inducing 
function of our identified miRNAs, we first identified the 
genes that have been known to relate to neuroblastoma 
differentiation using IPA. From this gene list, we identified 
the predicted targets for each miRNA and generated the 
predicted miRNAitarget interaction network mediating 
the differentiation-inducing function of the identified 
miRNAs. 

mRNA and miRNA expression. 

Total RNA was isolated as previously described 
[46]. mRNA expression profiling was done using the 
Illumina mRNA WG-6 v3 microarray platform. miRNA 
expression was measured by qRT-PCR using TaqMan 
microRNA Assays (Life Technologies) with average 
expression of RNU44 RNA, RNU19 RNA, GAPDH 
mRNA and 18s rRNA used as controls for normalizing 
RNA loading. 

Luciferase reporter assay. 

The segments of the wildtype 3 'UTRs for CDK4 and 
STATS containing the predicted target sites of miR-506-3p 
and miR-124-3p were cloned from human genomic DNA. 
Mutant constructs were generated with the seed target 
site (GUGCCUU) replaced by CACGGUU. The 3'UTRs 
were cloned downstream of the firefly luciferase coding 
sequences into the pmirGLO dual-luciferase reporter 
(Promega), a vector containing both firefly and Renilla 
luciferase cDNAs under the control of separate promoter/ 
terminator systems. The firefly luciferase was used as the 
primary reporter for miRNA regulation of the 3'UTR. The 
Renilla luciferase is an internal control for normalization. 
BE(2)-C cells were co-transfected with luciferase reporters 
(0.8 ng/ul) and miRNA mimics or control oligonucleotide 
(oligo) (25 nM). Luciferase activities were measured 
after 72 h using the Dual-Glo Luciferase Assay System 



(Promega). Firefly luciferase activity was normalized to 
Renilla luciferase activity to evaluate the effect of the 
miRNAs. 

Statistical analysis. 

For HCS, />-value for neurite length associated 
with each miRNA mimic was evaluated by assuming 
the distribution of neurite lengths follows a Gamma 
distribution, and FDR was determined by Benjamini- 
Hochberg correction method for multiple tests [47]. 
We consider a miRNA with FDR<0.01 as significantly 
inducing neurite outgrowth. For all other conditions, the 
statistical significance for each treatment was determined 
using ^test by comparing the treatment group with control, 
with /?< 0.05 considered statistically significant. 

ACKNOWLEDGEMENT 

The authors thank Dr. Luzhe Sun for critical reading 
of the manuscript. 

GRANT SUPPORT 

This work was supported in part by PRMRP 
Discovery Award PR121532 from the Department of 
Defense (to L. Du), NIH Cancer Center Support Grants 
P30 CA054174-17 (to L. Du, Y. Chen, A. Pertsemlidis, 
G.E. Tomlinson), NIH CTSA 1UL1RR025767-01 (to Y. 
Chen), and Greehey Children's Cancer Research Institute 
microRNA and Pediatric Cancer Pilot Award (to L. Du). 

REFERENCES 

1. Maris JM, Hogarty MD, Bagatell R, Cohn SL. 
Neuroblastoma. Lancet 2007;369(9579):2106-2120. 

2. Park JR, Eggert A, Caron H. Neuroblastoma: biology, 
prognosis, and treatment. Hematol Oncol Clin North Am 
2010;24(l):65-86. 

3. Brodeur GM. Neuroblastoma: biological insights into a 
clinical enigma. Nat Rev Cancer 2003;3(3):203-216. 

4. Reynolds CP. Differentiating agents in pediatric 
malignancies: retinoids in neuroblastoma. Curr Oncol Rep 
2000;2(6):51 1-518. 

5. Cruz FD, Matushansky I. Solid tumor differentiation 
therapy - is it possible? Oncotarget 2012;3(5):559-567. 

6. Nowak D, Stewart D, Koeffler HP. Differentiation 
therapy of leukemia: 3 decades of development. Blood 
2009;113(16):3655-3665. 

7. Matthay KK, Reynolds CP, Seeger RC, Shimada H, 
Adkins ES, Haas-Kogan D, Gerbing RB, London WB, 
Villablanca JG. Long-term results for children with 
high-risk neuroblastoma treated on a randomized trial 
of myeloablative therapy followed by 13-cis-retinoic 



www.impactjournals.com/oncotarget 



2510 



Oncotarget 



acid: a children's oncology group study. J Clin Oncol 
2009;27(7):1007-1013. 

8. Matthay KK, Villablanca JG, Seeger RC, Stram DO, Harris 
RE, Ramsay NK, Swift P, Shimada H, Black CT, Brodeur 
GM, Garbing RB, Reynolds CP. Treatment of high-risk 
neuroblastoma with intensive chemotherapy, radiotherapy, 
autologous bone marrow transplantation, and 13-cis- 
retinoic acid. Children's Cancer Group. N Engl J Med 
1999;341(16):1165-1173. 

9. Conrad C, Gerlich DW. Automated microscopy for high- 
content RNAi screening. J Cell Biol 2010;188(4):453-461. 

10. Jan E, Byrne SJ, Cuddihy M, Davies AM, Volkov Y, 
Gun'ko YK, Kotov NA. High-content screening as 
a universal tool for fingerprinting of cytotoxicity of 
nanoparticles. ACS Nano 2008;2(5):928-938. 

1 1 . Thiele C J, Reynolds CP, Israel MA. Decreased expression 
of N-myc precedes retinoic acid-induced morphological 
differentiation of human neuroblastoma. Nature 
1985;313(6001):404-406. 

12. Pahlman S, Odelstad L, Larsson E, Grotte G, Nilsson K. 
Phenotypic changes of human neuroblastoma cells in 
culture induced by 12-0-tetradecanoyl-phorbol- 13 -acetate. 
Int J Cancer 1981;28(5):583-589. 

13. Reynolds CP, Perez-Polo JR. Induction of neurite outgrowth 
in the IMR-32 human neuroblastoma cell line by nerve 
growth factor. J Neurosci Res 1981;6(3):319-325. 

14. Prasad KN, Kumar S. Role of cyclic AMP in differentiation 
of human neuroblastoma cells in culture. Cancer 
1975;36(4):1338-1343. 

15. Shenouda SK, Alahari SK. MicroRNA fiinction in cancer: 
oncogene or a tumor suppressor? Cancer Metastasis Rev 
2009;28(3-4):369-378. 

16. Li M, Li J, Ding X, He M, Cheng SY. microRNA and 
cancer. Aaps J 2010;12(3):309-317. 

17. Kinoshita T, Hanazawa T, Nohata N, Kikkawa N, Enokida 
H, Yoshino H, Yamasaki T, Hidaka H, Nakagawa M, 
Okamoto Y, Seki N. Tumor suppressive microRNA-218 
inhibits cancer cell migration and invasion through targeting 
laminin-332 in head and neck squamous cell carcinoma. 
Oncotarget 20 1 2;3( 1 1 ): 1 3 86- 1 400. 

18. Kopp F WE, Roidl A. The proto-oncogene KRAS is 
targeted by miR-200c. Oncotarget 2013. 

19. Bier A, Giladi N, Kronfeld N, Lee HK, Cazacu S, Finniss S, 
Xiang C, Poisson L, de Carvalho AC, Slavin S, Jacoby E, 
Yalon M, Toren A, et al. MicroRNA- 137 is downregulated 
in glioblastoma and inhibits the stemness of glioma stem 
cells by targeting RTVP-1. Oncotarget 2013;4(5):665-676. 

20. Kota J, Chivukula RR, O'Donnell KA, Wentzel EA, 
Montgomery CL, Hwang HW, Chang TC, Vivekanandan 
P, Torbenson M, Clark KR, Mendell JR, Mendell JT. 
Therapeutic microRNA delivery suppresses tumorigenesis 
in a murine liver cancer model. Cell 2009; 13 7(6): 1005- 
1017. 

2 1 . Kmtzfeldt J, Rajewsky N, Braich R, Rajeev KG, Tuschl T, 



Manoharan M, Stoffel M. Silencing of microRNAs in vivo 
with 'antagomirs'. Nature 2005;438(7068):685-689. 

22. Trang P, Wiggins JF, Daige CL, Cho C, Omotola M, Brown 
D, Weidhaas JB, Bader AG, Slack FJ. Systemic delivery 
of tumor suppressor microRNA mimics using a neutral 
lipid emulsion inhibits lung tumors in mice. Mol Ther 
2011;19(6):1116-1122. 

23. Lin RJ, Lin YC, Chen J, Kuo HH, Chen YY, Diccianni 
MB, London WB, Chang CH, Yu AL. microRNA signature 
and expression of Dicer and Drosha can predict prognosis 
and delineate risk groups in neuroblastoma. Cancer Res 
2010;70(20):7841-7850. 

24. Chang TC, Yu D, Lee YS, Wentzel EA, Arking DE, West 
KM, Dang CV, Thomas-Tikhonenko A, Mendell JT. 
Widespread microRNA repression by Myc contributes to 
tumorigenesis. Nat Genet 2008;40(l):43-50. 

25. Wei JS, Song YK, Durinck S, Chen QR, Cheuk AT, Tsang 
P, Zhang Q, Thiele CJ, Slack A, Shohet J, Khan J. The 
MYCN oncogene is a direct target of miR-34a. Oncogene 
2008;27(39):5204-5213. 

26. Tivnan A, Tracey L, Buckley PG, Alcock LC, Davidoff 
AM, Stallings RL. MicroRNA-34a is a potent tumor 
suppressor molecule in vivo in neuroblastoma. BMC 
Cancer 201 1;1 1:33. 

27. Makeyev EV, Zhang J, Carrasco MA, Maniatis T. The 
MicroRNA miR-124 promotes neuronal differentiation by 
triggering brain-specific alternative pre-mRNA splicing. 
Mol Cell 2007;27(3):435-448. 

28. Annibali D, Gioia U, Savino M, Laneve P, Caffarelli E, 
Nasi S. A new module in neural differentiation control: 
two microRNAs upregulated by retinoic acid, miR-9 and 
-103, target the differentiation inhibitor ID2. PLoS One 
2012;7(7):e40269. 

29. Le MT, Xie H, Zhou B, Chia PH, Rizk P, Um M, Udolph 
G, Yang H, Lim B, Lodish HF. MicroRNA- 125b promotes 
neuronal differentiation in human cells by repressing 
multiple targets. Mol Cell Biol 2009;29(19):5290-5305. 

30. Foley NH, Bray I, Watters KM, Das S, Bryan K, Bernas T, 
Prehn JH, Stallings RL. MicroRNAs 10a and 10b are potent 
inducers of neuroblastoma cell differentiation through 
targeting of nuclear receptor corepressor 2. Cell Death 
Differ 2011;18(7):1089-1098. 

3 1 . Swarbrick A, Woods SL, Shaw A, Balakrishnan A, Phua 
Y, Nguyen A, Chanthery Y, Lim L, Ashton LJ, Judson 
RL, Huskey N, Blelloch R, Haber M, et al. miR-380-5p 
represses p53 to control cellular survival and is associated 
with poor outcome in MYCN-amplified neuroblastoma. Nat 
Med 2010;16(10):1134-1140. 

32. Mishra PJ, Merlino G. MicroRNA reexpression 
as differentiation therapy in cancer. J Clin Invest 
2009;119(8):21 19-2123. 

33. Radio NM, Breier JM, Shafer TJ, Mundy WR. Assessment 
of chemical effects on neurite outgrowth in PC 12 cells 
using high content screening. Toxicol Sci 2008; 105(1 ):1 06- 



www.impactjournals.com/oncotarget 



2511 



Oncotarget 



118. 

34. Cheung YT, Lau WK, Yu MS, Lai CS, Yeung SC, So KF, 
Chang RC. Effects of all-trans-retinoic acid on human SH- 
SY5Y neuroblastoma as in vitro model in neurotoxicity 
research. Neurotoxicology 2009;30(1): 127-135. 

35. Mao L, Ding J, Zha Y, Yang L, McCarthy BA, King W, 
Cui H, Ding HF. H0XC9 Unks cell-cycle exit and neuronal 
differentiation and is a prognostic marker in neuroblastoma. 
Cancer Res 20 1 1 ;7 1 ( 1 2):43 1 4-4324. 

36. Smith B, Treadwell J, Zhang D, Ly D, McKinnell I, Walker 
PR, Sikorska M. Large-scale expression analysis reveals 
distinct microRNA profiles at different stages of human 
neurodevelopment. PLoS One 2010;5(6):el 1 109. 

37. Silber J, Lim DA, Petritsch C, Persson AT, Maunakea AK, 
Yu M, Vandenberg SR, Ginzinger DG, James CD, Costello 
JF, Bergers G, Weiss WA, Alvarez-Buylla A, et al. miR- 
124 and miR-137 inhibit proliferation of glioblastoma 
multiforme cells and induce differentiation of brain tumor 
stem cells. BMC Med 2008;6:14. 

38. Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein 
EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, 
Rajewsky N. Combinatorial microRNA target predictions. 
Nat Genet 2005;37(5):495-500. 

39. Lewis BP, Burge CB, Bartel DP. Conserved seed pairing, 
often flanked by adenosines, indicates that thousands of 
human genes are microRNA targets. Cell 2005 ; 120(1) : 15- 
20. 

40. Stark A, Brennecke J, Bushati N, RusseU RB, Cohen SM. 
Animal MicroRNAs confer robustness to gene expression 
and have a significant impact on 3'UTR evolution. Cell 
2005;123(6):1 133-1 146. 

41. Laurent LC, Chen J, Ulitsky I, MueUer FJ, Lu C, Shamir 
R, Fan JB, Loring JF. Comprehensive microRNA 
profiling reveals a unique human embryonic stem cell 
signature dominated by a single seed sequence. Stem Cells 
2008;26(6):1506-1516. 

42. Yeyeodu ST, Witherspoon SM, Gilyazova N, Ibeanu GC. 
A rapid, inexpensive high throughput screen method for 
neurite outgrowth. Curr Chem Genomics 2010;4:74-83. 

43. Price RD, Oe T, Yamaji T, Matsuoka N. A simple, flexible, 
nonfluorescent system for the automated screening of 
neurite outgrowth. J Biomol Screen 2006; 1 1(2): 155-164. 

44. Mitchell PJ, Hanson JC, Quets-Nguyen AT, Bergeron M, 
Smith RC. A quantitative method for analysis of in vitro 
neurite outgrowth. J Neurosci Methods 2007;164(2):350- 
362. 

45. Liang Y, Ridzon D, Wong L, Chen C. Characterization of 
microRNA expression profiles in normal human tissues. 
BMC Genomics 2007;8:166. 

46. Du L, Schageman JJ, Subauste MC, Saber B, Hammond 
SM, Prudkin L, Wistuba, II, Ji L, Roth JA, Minna JD, 
Pertsemlidis A. miR-93, miR-98, and miR-197 regulate 
expression of tumor suppressor gene FUSl. Mol Cancer 
Res 2009;7(8): 1234-1243. 



47. Yoav Benjamini YH. Controlling the False Discovery 
rate: A Practical and Powerful Approach to Multiple 
Testing. Journal of the Royal Statistical Society Series B 
(Methodological) 1995;57(1):12. 



www.impactjournals.com/oncotarget 



2512 



Oncotarget 



